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INTRODUCTION 
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1.1  Rationale  for  the  Measurements/Calculations 

This  document  summarizes  the  neutron  spectral  measurements  and  the 
neutron  and  gamma-ray  spectral  calculations  for  five  different  irradiation  geo¬ 
metries  commonly  employed  in  biological  and  semiconductor  radiation-damage  studies 
free  field  in  both  ER1  and  ER2  exposure  rooms;  behind  a  2-inch  lead  wall  and 
inside  an  exercise  wheel;  behind  a  6-inch  lead  wall  and  inside  a  2-inch  lead 
cave;  and  at  the  center  of  a  tissue-equivalent  phantom  located  behind  a  6- 
inch  lead  wall.  Not  only  does  the  gamma  to  neutron  dose  ratio  change  consid¬ 
erably  from  one  of  these  configurations  to  the  next,  but  the  neutron  spectrum 
changes  very  significantly  for  each  perturbation  of  the  free-field  measurement. 

(In  fact,  the  free-field  spectrum  itself  changes  measurably  with  distance 
from  the  reactor1,  because  of  the  relatively  greater  effect  of  room-scattered 
neutrons  with  increasing  distance  (R)  from  the  center  of  the  reactor:  The 
fast  neutrons  from  the  reactor  fall-off  approximately  as  1/R2  while  the  room- 
scattered  component  Is  more  nearly  constant  with  R,  the  distance  from  the  re¬ 
actor  core  center.)  The  methods  of  "hardening"  or  "softening"  the  fast  neu¬ 
tron  spectrum  are  presented  In  Section  4. 

Accurate  neutron  spectral  measurements  were  made  for  the  five  dif¬ 
ferent  configurations.  The  neutron-transport  calculations  were  carried  out 
for  the  same  cases  with  geometric  approximations  that  varied  from  a  spherical 
geometry  (of  various  "shells"  of  reactor  core,  aluminum,  air,  lead,  water, 
and  the  distant  walls)  to  a  3-dimensional  (3-D)  model  of  the  actual  geometry. 

The  measurements  utilized  the  calculated  spectrum  shape  in  each  case  as  a 
starting  (trial)  spectrum,  and  the  final  unfolded  spectrum  was  obtained  with 
only  a  few  small  shape-varying  Iterations.  This  only  occurs  when  the  cal¬ 
culated  spectrum  and  the  measured  are  in  good  agreement,  and  verifies  the 
accuracy  of  the  calculated  spectral  shape. 
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Once  the  accuracy  of  the  calculated  spectral  shapes  have  been 
verified  with  multiple-foil  measurements,  the  magnitude  Is  obtained  with 
the  use  of  a  single  multiplying  factor  (normalization  constant).  This 
normalization  factor  is  required  to  correct  for  uncertainties  in  the  abso¬ 
lute  neutron  output  of  the  reactor  per  power-monitor  count  (power  calibration), 
in  power  distribution  throughout  the  reactor,  and  the  (variable)  water  thick¬ 
ness  layer  between  the  reactor  core  and  the  exposure-room  aluminum  "window". 
Some  smaller  uncertainties  can  arise  from  cross-section  approximations  and 
finite  energy  group  widths  but  these  and  the  previously  mentioned  uncertainties 
are  completely  resolved  by  renormalization  of  the  calculation  (by  determining 
the  simple  "correction  constant")  to  the  nearest  configuration  measured. 

With  the  baseline  measurements/calculations  presented  here,  one 
can  very  accurately  calculate  the  biological  dose  and  thus  calibrate  all  the 
fast  neutron  dosimetry  devices  used  with  the  AFRRI  reactor,  inluding  the 
sulfur  or  nickel  neutron  flux  monitor.  All  of  the  above  applies  to  fast 
neutron  flux. 

For  thermal  neutrons,  a  cadmium-difference  measurement  with  manga¬ 
nese  or  (very  thin)  gold  foils  or  other  thermal -neutron  sensitive  device  is 
adequate.  The  measurement  of  the  thermal  neutron  flux  is  an  absolute,  single- 
point  measurement.  On  the  other  hand,  the  accurate  calculation  of  thermal- 
neutron  flux  is  subject  to  variations  due  to  small  material  and  geometric  un¬ 
certainties,  including  the  amount  of  neutron-absorbing  impurities  in  the  room 
structure  and  that  of  the  added  shields/moderators  which  is  generally  very 
poorly  known.  Thus, a  careful  measurement  is  much  preferred. 

Returning  to  the  fast  neutron  flux,  it  should  be  pointed  out  that 
once  the  spectral  shape  has  been  accurately  determined  by  a  calculation  such 
as  the  ones  experimentally  verified  here,  only  one  threshold  foil  dosimeter 
is  subsequently  needed  to  fix  the  absolute  magnitude:  Thus  the  power  of  .the 
calculations  for  future  use  such  as  in  the  design  of  biological/material 
exposure  is  clear. 
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The  gamma-ray  dose  variability  can  be  ascertained  with  adequate 
accuracy  using  a  number  of  well  characterized  dosimeters.  The  RBE  (relative 
biological  effectiveness)  does  not  depart  much  from  unity  for  gamma-rays, 
unlike  neutrons  for  which  the  RBE  varies  by  about  a  factor  of  10  between 
slow  and  fast  neutrons.  The  neutron  damage  to  electronic  devices  varies  by 
even  a  greater  factor2.  Thus,  for  the  AFRRI  reactor  exposures,  the  charac¬ 
terization  of  the  radiation  field  in  the  work  presented  below  is  concerned 
chiefly  with  determinations  of  the  neutron  spectral  flux. 

The  calculations  of  gamna-ray  spectra  (presented  in  Volume  II  of 
this  report)  are  obtained  by  using  coupled  neutron  and  gamma-ray  cross  sec¬ 
tions.  The  gamma-ray  field  is  mostly  made  up  of  direct  gamma -rays  from  the 
reactor  core,  inelastic  scattering  of  fast  neutrons,  and  thermal -neutron- 
capture  gamma-rays.  Because  of  the  uncertainty  in  thermal  neutron  flux,  the 
calculation  of  the  capture  gamma-ray  flux  is  not  expected  to  be  nearly  as 
accurate  as  the  calculation  of  the  fast  neutron  flux:  This  is  important  be¬ 
hind  lead  shields  where  a  large  fraction  of  the  total  gamna-ray .dose  may  be 
due  to  capture-neutron  ganma-rays. 

1.2  Past  Work 

The  spectral  characteristics  of  the  neutron  field  were  measured 
almost  two  decades  ago  by  W.  Quam  of  E.G.&G,  and  first  reported  in  October, 
1964.  A  later  summary  appeared  in  a  compilation  by  J.A.  Sayeg3.  In  this  work, 
the  gross  spectral  dependence  of  the  neutron  field  was  measured  with  the  use 
of  three  threshold  foils,  and  also  with  some  foils  without  an  energy  threshold 
(Eth)  but  with  boron  and  cadmium  shields  to  produce  effective  (mechanical)  low 
energy  thresholds  of  about  0.01  MeV  and  4xl0“7  MeV,  respectively.  The  foils 
and  their  threshold  energies  were  32S(Eth  *  3  MeV),  2S8U(Eth  *  1.5  MeV), 
237Np(E^  0*5  MeV),  2,*Pu(EtjJ  ^  0-01  MeV  via  boron  shield),  and  1,7Au(E^ 

0.01  MeV,  4xl0"7  MeV,  and  0  MeV  via  boron  shield,  cadmium  shield,  and  un¬ 
shielded  foils  respectively.) 

For  the  fission  foils,  gross  fission-product  gamma-ray  counting  was 
done,  with  essentially  all  the  gamma  rays  being  counted.  The  count  rate  versus 
irradiation  time,  wait  time,  and  count  time  was  calibrated  with  a  very  complex 
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procedure  which  leaves  much  to  be  desired  in  terms  of  accuracy.  Furthermore, 
the  non-fission- product  radioactivity  must  be  somehow  estimated  and  subtracted 
in  this  type  of  measurement.  This  adds  significantly  to  the  uncertainties  of 
these  early  measurements,  especially  for  fission  foils  such  as  237Np  which  are 
highly  radioactive  before  being  exposed  to  the  neutron  field  being  measured. 

1.3  Present  Work 

The  state  of  the  art  has  advanced  considerably  since  the  early  1960's 
in  neutron  spectrum  measurements  in  terms  of  radioactivity-counting  techniques; 
nuclear  decay  schemes/fluorescent  yields/fission-product  yields;  cross  section 
data;  large-computer-based  spectrum-unfolding  codes;  and  neutron  transport  cal¬ 
culations  to  provide  a  good  trial  spectrum  for  the  neu ;ron-spectrum  unfolding 
code.  These  advancements  are  presented  in  recent  ASTM  Procedures'*,  which  were 
essentially  derived  from  the  work  presented  in  Reference  5.  The  work  presented 
here  was  carried  out  by  the  same  methods  and,  in  fact,  by  the  same  investigators. 
Some  recent  improvements  have  been  added.  These  involve  some  more  accurate 
nuclear  data  on  the  threshold-foil  cross  sections  and  gamma-ray  fields,  and 
more  careful  measurements  of  thermal  and  near-epithermal  neutron  spectral  data 
(which  were  of  little  interest  to  radiation  damage  characterization  to  silicon 
semiconductor  devices  addressed  in  the  ASTM  Standard  Procedures4). 
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2. 


FIELD  CHARACTERIZATION  MEASUREMENTS 


2. 1  Foil  Measurement  Methods 

Five  baseline  spectral  measurements  were  performed  corresponding  to 
the  most  frequently  utilized  irradiation  configurations.  The  characteristics 
of  these  five  reactor  runs  are  shown  in  Table  1  and  the  experimental  confi¬ 
guration  for  each  is  shown  in  Figures  1  through  5.  The  characteristics  of  the 
seventeen  threshold  foils  used  to  measure  the  above  neutron  spectra  are  shown 
in  Table  2. 

Three  threshold  foil  sets  were  used  for  each  spectral  measurement. 

The  first  was  a  complete  set  of  all  seventeen  foils  and  was  enclosed 
in  a  1.65  g/cm2  boron  sphere;  the  second  set  contained  one  1,7Au,  one  #5Mn  and 
three  S8Ni  foils,  all  enclosed  in  a  0.020  inch  (.050  cm)  thick  cadmium  cover, 
and  the  third  set  contained  one  197Au,  one  S5Mn  and  three  5#Ni  foils,  with  all 
five  foils  in  "bare"  (no  cover)  configuration. 

The  use  of  the  ,0B  sphere  is  important  for  two  reasons.  First,  it 
provides  an  effective  threshold  of  about  0.01  MeV  for  the  197Au,  the  55Mn, 
and  2J5U  foils;  and  second,  it  reduces  drastically  the  effects  of  isotopic 
impurities  in  the  foils  that  lead  to  the  same  product  nucleus  as  the  reactions 
of  interest  via  the  much  more  efficient  route  of  thermal  neutron  activation. 

For  instance,  the  reactions  Al(n,  a)  and  2SNa(n,  y)  lead  to  the  same  product 
nucleus,  but  the  second  one  is  activated  by  thermal  neutrons  that  have  a  very 
large  reaction  cross  section.  Thus,  in  this  case,  sodium  contamination  of  the 
aluminum  foils  (surface  contamination,  say)  must  be  kept  very  small. 

The  use  of  the  cadmium-covered  and  "bare"  foils  provides  an  excellent 
method  for  determining  the  thermal  neutron  flux  by  first  determining  the  total 
flux  from  the  "bare"  foils  and  then  subtracting  out  the  epithermal  flux  determined 
from  the  cadmium-covered  foils. 

Finally,  the  58N1  foils  in  the  "bare"  foil  set  are  used  in  conjunction 
with  the  **Ni  foils  in  the  1#B  sphere  to  provide  a  correction  factor  for  the 
,08  shielding  effect  on  the  foils  with  Eth  £  1  MeV.  These  s,Ni  foils  are  used  with 
*2S  foils  at  some  dosimetry  laboratories  because  they  are  Inherently  more  accurate. 
This  stems  from  utilizing  a  monoenergetic  gamma  ray,  versus  beta  counting  for  ,2S 
activation  measurement. 
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Table  1.  AFRRI-Reactor-Run  Characteristics. 


Run  # 

Exp.  Room 

Run  Designation  Power  Level 

(kW) 

Run  Duration 
(Sec. ) 

1 

ER-1 

6"  Lead  Wall  /  Cave 

250 

1800 

2 

ER-2 

Free-Field 

250 

600 

3 

ER-1 

Exercise  Wheel 

250 

900 

4 

ER-1 

Free-Field 

250 

600 

5 

ER-1 

Phantom  /  6"  Lead  Wall 

500 

1200 

! 

t 
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SAI-8 1CC-20.  1 


Figure  1.  Top-view  Schematic  Diagram  of  AFRRI 
Reactor  Run  #1  (ER1  -  6"  Lead  Wall  +  2” 
Lead  Cave). 


SAI-81CC-2J 


Figure  2.  Schematic  Diagram  of  AFRRI 
Reactor  Run  12  {ER-2  Free  Field).  Top  View. 
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REACTOR  CORE  CENTER 


ER-1 


SAI-81cc-24 

Flqure  3.  Schematic  Diagram  of  AFRRI 
Reactor  Run  #3  (ER-1  EXERCISE  WHEEL  behind  2"  Lead  Wall).  Top  View. 


Figure  4.  Schematic  Diagram  of  AFRRI 
Reactor  Run  #4  (ER-1  Free  Field).  Top  View. 


i  « 
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SAI-81CC-23. 1 


Figure  5.  Top-view  Schematic  Diagram  of  AFRRI 
Reactor  Run  *5  (ER-1  "PHANTOM"  behind  6"  Lead  Wall). 


The  AFRRI  Ge(Li)  detector  was  used  for  foil  activation  counting  in 
two  geometric  configurations.  The  hotter  foils  were  placed  farther  from  the 
detector  to  eliminate  excessive  dead-time  effects.  The  detector  was  calibrated 
for  both  geometries  using  an  NBS  eleven-line  mixed  radionuclide  source.  Some  of 
the  longer  lived  foils  were  counted  at  the  SAI  Ge (Li)  detector. 

The  live-time  for  each  foil  count  was  accurately  determined  by  using  a 
calibrated  pulser  as  the  test  input  to  the  detector  preamplifier:  Any  pulser 
event  that  occurred  during  the  time  the  analyzer  was  busy  or  during  incomplete 
recovery  of  the  amplifier  baseline  from  a  previous  gamma-ray  pulse  would  cause 
the  pulser  pulse  to  lose  a  count  in  the  (high  resolution)  peak. 

At  the  conclusion  of  all  foil  counting  the  most  active  foil 
was  used  to  determine  the  cascade  correction  to  be  applied  on  the  fission 
foil  activations  as  determined  from  the  537  and  1596  keV  lines.  The  cascade 
effect  for  these  particular  ganma-rays  is  fully  discussed  in  Ref.  1  and  it 
affects  the  detector  efficiency  at  these  energies.  Its  determination  consists 
of  measuring  the  relevant  count  rates  R^  at  different  distances  from  the  detector 
face  and  plotting  as  a  function  of  distance  d  from  the  detector  housing. 

The  fractional  difference  of  1//r. at  d  *  0  (the  surface  of  the  detector 
cryostat),  from  the  straight  line  plot  for  d  >  0  values,  gives  a  measure  of 
the  magnitude  of  the  cascade  correction  for  the  particular  detector  used.  It 
was  found  that  this  correction  was  of  the  order  of  15%  for  the  AFRRI  Ge(Li) 
detector  at  537  and  1596  keV. 

The  activations  per  nucleus  induced  in  each  foil  were  determined 
from  the  measured  foil  counts  by  the  following  equations 


N  =  NA  Xti  eXt»'(l-e'Xti)'1(l-e~uc)'1 
eYfYyN  m(a/A) 


(1) 


where 


is  the  area  of  the  gamma-ray  peak 

X  is  the  decay  constant  of  the  reaction 

tj  is  the  irradiation  time 

Is  the  time  interval  between  the  end  of  irradiation 
"  and  the  beginning  of  counting 

tc  Is  the  foil  count- time 

c  is  the  detector  efficiency 

Yy  is  the  gamna-ray  yield  per  reaction 


Yf  is  the  fractional  fission  yield  (set  this  to  unity 
No  is  Avogadro's  number  for  non-fission  foils), 

m  is  the  foil  mass 

a  is  the  fractional  natural  abundance  of  the  reactive 
isotope  under  consideration,  and 

A  is  the  atomic  weight  of  the  foil  element. 


140  140 

For  the  1596  keV  line  of  La  which  is  produced  by  the  Ba  decays. 
Equation  (1)  becomes 

N  =  Na  XBti(XL-XB)  (l-e-*8^)'1  (i-e~ABtc)“1  (2) 

>LeY^f(e"XBtw-e'XLtw)Nom(a/A) 

where  the  subscripts  (superscripts)  B  and  L  refer  to  quantities  associated 
140  140 

with  the  Ba  and  La  garanas  respectively. 

Finally,  as  a  cursory  check  against  possible  errors  in  determining 
the  number  or  activations  per  nucleus  from  Equations  (1)  and  (2)  above,  the 
ratio  N/Nn-  is  determined  for  all  foils  in  the  reactor  run  under  consideration, 
and  compared  to  N/N^  from  past  experiments  on  file  with  similarly  configured 
reactors.  (Here,  Nw.  is  the  number  of  activations  per  nucleus  induced  in  the 
nickel  foil  by  the  reaction  (  Ni(n,  p)30Co).)  In  addition,  this  excercise 
provides  first-hand  qualitative  information  on  the  approximate  shape  of  the 
neutron  spectrum,  prior  to  the  final  unfolding  process  which  is  discussed 
in  the  next  section. 

The  foil  activation  data  determined  from  reactor  Runs  1  through  5  are 
presented  in  Table  2,  and  the  unfolded  spectra  are  shown  in  Figures  6 
through  10. 
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Table  2.  THRESHOLD  FOIL  ACTIVATIONS**  FOR  AFRRI  TRIGA 
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Aborted  readings  (inconsistent  with  the  other  foils). 

Read  as  2.976  x  10"®. 

The  l^Ba,  ^Zr,  an<j  140|_a  lines  are  used  in  conjunction  with  the 
fission  foils,  except  for  232Th,  where  !4°la  cannot  be  used. 
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Figure  8.  E*(j>f  (E)  for  AFRRI  Reactor  Run  #3  (Exercise  Wheel  behind  2"  Pb  Wall). 


(E)  for  AFRRI  Reactor  #4  (ER-1 


Figure  10.  E*<K  (E)  for  AFRRI  Reactor  #5  (Phantom  behind  6"  Pb  Wall). 


2.2  Results:  Shape-Comparison  Between  Measured 

and  Calculated  (Trial)  Spectra 

The  unfolded  neutron  spectra  resulting  from  the  measured  activations 
are  shown  in  Figures  11  through  15  for  reactor  Runs  1  through  5,  respectively. 
Also  shown  are  the  calculated  spectra  for  the  same  reactor  configurations, 
which  were  used  as  the  input  trial  for  the  SAND  II  code.  For  shape  comparison 
the  measured  and  calculated  spectra  are  normalized  so  that  best  fit  is  obtained 
for  the  energy  region  above  En  ^0.2  MeV,  which  is  the  effective  threshold  for 
semiconductor  neutron  damage.  The  measured  spectra  are  normalized  to  unit 
fluence  above  10  keV. 

With  regard  to  the  calculated  trial  spectra  it  must  be  noted  that 
both  one-dimensional  and  three-dimensional  transport  codes  were  utilized  in 
modeling  the  actual  experimental  conditions.  In  Runs  1,  2  and  4  the  one-di¬ 
mensional  models  were  found  to  be  adequate  as  trial  spectra  in  the  unfolding 
code. 

For  reactor  Run  #3  (exercise  wheel  behind  2"  Pb  shield)  and  Run  #5 
(phantom  behind  6"  Pb  shield),  it  was  found  that  a  mixed  1-D/3-D  trial  spectrum 
produced  much  better  shape  agreement  with  the  measured  ones  than  either  the 
1-D  or  3-D  spectrum  alone.  The  "mixing"  was  accomplished  by  splicing  the 
high  energy  part  (E^4MeV)  of  the  1-D  calculation  to  the  3-D  calculation  part 
with  E<4MeV,  where  the  latter  does  not  suffer  much  due  to  the  statistical 
uncertainity.  The  need  for  a  spliced  spectrum  trial  spectrum  in  Run  3  (2"  Pb 
wall  and  activity  wheel)  and  Run  5  (phantom  behind  6"  Pb  shield)  appears  due 
to  a  combination  of  reasons.  The  complex  geometry  of  the  activity  wheel  and 
the  large  thermalizing  power  and  capture  properties  of  the  lucite  wheel  and 
tissue-equivalent  phantom  required  the  need  for  the  3-D  calculations.  Un¬ 
fortunately,  as  stated  above,  the  statistical  uncertainty  of  the  high  energy 
groups  in  the  3-D  spectra  required  replacement  by  the  1-D  spectra.  Volume  II 
of  this  report  contains  further  discussion  of  the  3-D  and  spliced  calculations. 

It  is  seen  that  the  calculated  and  the  measured  spectra  for  reactor 
Runs  1  through  4  agree  quite  well  in  shape:  It  took  no  more  than  two  iterations 
for  the  unfolding  code  to  reach  the  required  5%  deviation  between  measured 
and  calculated  foil  activations. 
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The  measured  data  were  normalised  to  unit  flux  above  20  keV  and  the 
calculated  data  normalized  to  the  measured  data  above  0.2  meV,  the 
"threshold"  for  radiation  damaye  to  silicon. 


In  reactor  Run  5  (Phantom),  the  measured  and  calculated  spectra  do  not  agree 
well  in  shape  in  that  the  calculation  yields  about  a  50%  excess  in  the  1/E 
region  of  the  spectrum  (below  0.3  MeV).  As  mentioned  above,  the  calculated 
spectrum  for  this  case  is  a  mixture  of  a  one-dimensional  and  a  three-dimensional 
calculation  joined  at  En  =  4.07  MeV.  In  the  neighborhood  of  this  energy,  the 
3-D  calculation  still  suffers  from  statistical  uncertainties.  Thus,  the  lack 
of  shape  agreement  in  this  case  is  at  least  partly  due  to  difficulties  in 
accurate  normalization  between  the  two  calculations  at  their  juncture  point. 

The  tabulated  spectra  for  Runs  1  through  5  are  presented  in  Table  3. 
These  data  represent  the  absolute  flux,  per  kilowatt,  for  the  calculated  data 
as  well  as  for  the  measurements:  There  is  no  normalization  to  the  measured 
data  above  0.2  MeV,  as  in  Figures  11  through  15,  so  that  direct  comparisons 
can  be  made  as  to  predicted  group  flux  per  kilowatt. 
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BROAD-GROUP  COMPARISON  BETWEEN  MEASUREMENTS  AND  CALCULATIONS 


The  measured  neutron  spectra  shown  in  Figures  11  thourgh  15  were 
regrouped  in  the  same  37  energy-group  structure  as  the  calculations  for 
point-by-point  comparison.  The  results  are  presented  in  Table  3.  The  ab¬ 
solute  measured  and  calculated  data  are  shown;  i.e.,  no  normalization  was 
employed  here.  The  calculations  can  most  meaningfully  be  compared  to  the 
measured  data  in  the  0.1  to  7  MeV  energy  region,  where  most  of  the  biological 
dose/silicon  damage  occurs. 

Table  4  presents  measured  and  calculated  fluxes  above  10  keV, 
above  the  cadmium  cutoff  (0.4  eV),  and  for  the  thermal  group  (below  0.4  eV). 
Agreement  is  excellent  for  ER-2  free  field  and  is  good  for  the  6"  Pb  plus 
2"  cave.  In  all  cases  except  Run  5  (phantom  behind  6"  of  lead)  agreement 
between  the  measured  versus  calculated  fast  and  epicadmium  fluxes  (ratios 
between  measured  groups  1-27  and  1-36  and  the  corresponding  calculated  groups) 
is  quite  good.  In  a  similar  fashion  the  agreement  between  measured  and  cal¬ 
culated  thermal  fluxes  is  poor  when  compared  to  the  agreement  between  the 
fast  and  the  epicadmium- region  fluxes.  It  is  to  be  noted  that  the  spliced 
1-D/3-D  spectra,  (Runs  3  and  5)  have  the  greatest  differences.  These  dif¬ 
ferences  may,  in  part,  be  ascribed  to  a  renormalization  difference  between 
the  1-D  and  3-D  calculations.  In  general  the  calculation  of  thermal  neutron 
flux  is  far  more  complicated  than  for  fast  neutrons  because  thermal  neutrons 
are  very  sensitive  to  the  thickness  and  location  of  all  scattering  materials, 
moderating  materials  (wood,  plastic,  concrete  etc.),  and  of  the  thermal 
neutron  absorbers. 

4.0  SPECTRUM  TAILORING  WITH  THE  AFRRI  NEUTRON  FIELD 

In  comparing  the  measured  free  field  flux  (Run  #4,  ER-1)  and  the 
flux  behind  the  6"  lead  wall  with  2"  lead  cave  (Run  #1),  note  that  the  flux 
is  higher  behind  the  6"  lead  wall  with  2"  lead  cave  than  in  front  of  the  wall 
(Run  #4  free  field)  in  the  0.1  to  1.0  MeV  energy  region.  Above  1  MeV,  the 
flux  behind  the  wall  is  appreciably  lower.  A  comparison  of  the  spectral 
shapes  of  Figures  11  (6"  Pb  and  cave)  and  12  (free  field)  shows  the  reason 
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Table  4.  Absolute  Comparison  Between  Calculated  and  Measured  Flux  Per 
Kilowatt  for  Three  Spectral  Regions. 


Groups  1-27  Groups  1-36  Group  37 

Run  §  (En  >  10  keV)  (Epicadmium  Region)  (Thermal  Region) 


- , 

Hi 

Meas. 

H 

<1>  Meas. 

<t>  Meas. 

1 

6.85+7 

7.88+7 

9.22+7 

1.04+8 

2.79+6 

2.63+6 

2 

7.83+7 

7.89+7 

1.17+8 

1.22+8 

3.37+7 

5.59+7 

3 

1.06+8 

8.20+7 

1.47+8 

1.18+8 

5.55+7 

2.06+7 

4 

9.86+7 

7.66+7 

1.32+8 

1.08+8 

3. 34+6 

3.61+6 

5 

2.08+7 

1.91+7 

3.77+7 

3.02+7 

1.02+8 

5.58+7 

for  this:  The  6"  lead  wall  downscatters  the  high  energy  neutrons  by  a  sig¬ 
nificant  amount,  with  the  result  that  the  E*$(E)  "spectrum"  for  the  6"  Pb 
penetration  peaks  at  about  0.8  MeV,  whereas  the  corresDonding  source  "spec¬ 
trum"  (Figure  12)  peaks  at  about  3  MeV.  This  is  due  to  the  energy  dependence 
of  the  neutron  inelastic  scattering  cross  section  for  lead;  it  is  roughly 
equal  to  50%  of  the  total  cross  section  at  6  MeV  and  above,  and  drops  by  over 
an  order  of  magnitude  from  6  MeV  to  1  MeV,  vrtiere  it  becomes  essentially  negli¬ 
gible.  Thus,  the  neutrons  below  1  MeV  are  "built  up"  by  down-scattering  from 
higher  energies.  The  result  is  a  shifting  of  the  peak  of  the  E*$(E)  "spectrum" 
from  about  3  MeV  down  to  about  0.8  MeV.  (See  Figures  11  and  14  above.)  If 
one  normalizes  the  two  spectra  in  the  0.1  to  0.8  MeV  region,  the  two  will 
be  seen  to  differ  by  a  factor  of  3  to  4  in  the  few  MeV  region.  In  previous 
work  done  by  the  authors2,  this  degree  of  spectrum  "softening"  was  obtained 
with  only  2"  of  iron  shielding:  The  inelastic  scattering  cross  section  for 
iron  is  quite  large  down  to  about  0.8  MeV.  Thus,  with  a  medium  to  high  atomic 
weight  shield,  the  neutron  spectrum  becomes  appreciably  softened. 

It  is  interesting  to  note  that  the  spectral  difference  between  the 
free  field  and  6i!  Pb  wall  +  cave  cases  is  great  enough  wherein  the  thought  of 
useful  spectral -dependence  experiments  might  be  entertained  for  comparing 
"equidose"  exposures  for  the  two  cases  of  widely  differing  neutron  spectra. 

Such  experiments  would  probably  have  to  be  confined  to  those  where  the  gamma- 
ray  sensitivity  is  low  because  it  is  much  different  relative  to  the  two  neu¬ 
tron  dosages.  An  example  of  this  is  the  measurement  of  transistor  damage 
in  different  neutron  spectra. 

With  a  water  shield,  the  effect  is  a  hardening  of  the  neutron 
spectrum.  This  arises  from  a  rapid  overall  decrease  in  neutron  cross  section 
with  increasing  energy  for  both  hydrogen  and  oxygen  above  1  MeV:  The  higher 
energy  neutrons  become  quite  transparent.  A  comparison  of  the  free-field. 
spectrum.  Figure  16,  with  the  spectrum  calculated  for  the  reactor  moved  back 
12"  into  the  water  (Figure  17)  shows  a  dramatic  difference.  The  peak  energy 
for  the  12"  water  shield  is  moved  up  to  the  4-8  MeV  region.  Also,  the  E*$(E) 
"spectrum"  drops  by  about  two  orders  of  magnitude  below  the  peak  value  for 
energies  below  0.1  MeV  where  the  1/E  spectral  dependence  finally  sets  in. 
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(Front+Back)  Neutron  Flux  vs  Energy  ER1  Free  Field  with  12"  Water 


The  above  methods  of  obtaining  very  diverse  spectra  were  used  for 
2 

experiments  with  transistors  ,  wherein  the  measured  damage  was  compared  to  that 
calculated  for  exposure  to  widely  varying  neutron  spectra.  The  validity  of  the 
calculated  silicon  damage  function  D( En )  was  verified,  since  the  measured  damage 
and  that  calculated  from  ^D( En (En )dEn  agreed  for  exposures  to  widely  differing 
spectra.  With  the  use  of  only  2"  of  iron  to  obtain  a  softened  fast-neutron 
spectrum,  and  4"  of  polyethylene  for  a  hardened  one,  the  two  differed  by  about 
a  factor  of  10  above  several  MeV,  when  normalized  at  the  0.2  MeV  threshold  energy 
for  damage  to  silicon  semi  conductor  material. 

In  this  way,  a  TRIGA  reactor  such  as  the  AFRRI  reactor  can  be 
used  to  study  the  energy  dependence  of  other  semiconductor  materials  such 
as  GaAs  or  Ge,  of  electrical  insulators  materials,  of  structural  materials 
such  as  steel,  and  other  radiation  effects  which  do  not  have  a  strong  gamma- 
ray  irradiation  dependence. 


5.  CHARACTERIZATION  OF  THE  NEUTRON  FIELDS  FOR  SILICON  DAMAGE 

EFFECTIVENESS 

The  silicon-damage  characterization  for  each  of  the  fi.e  measure¬ 
ments  is  presented  in  Table  5.  These  values  are  for  total  fluence  and  not 
flux/kW.  Note  in  Table  1  that  the  free  field  measurements  (Runs  2  and  4) 
were  carried  out  fora  total  of  150  MW-sec,  the  6"  lead  wall  +  2"  cave  case 
(Run  1)  for  450  MW-sec,  the  exercise  wheel  (2"  Pb)  for  225  WM-sec,  and  the 
phantom  behind  6"  Pb  for  600  MW-sec. 

Normalizing  to  a  standard  factor  power  *  time,  one  finds  that  the  1 
MeV  equivalent  flux  (<t*eq»  Column  5)  is  28%  higher  for  the  free  field  case 
than  for  the  6"  Pb  wall  shield  and  2"  cave.  Yet,  the  flux  above  10  keV  is 
the  same  for  both.  Thus,  the  softer  neutrons  emerging  from  the  6"  lead  wall 
have  28%  less  effect  on  silicon.  The  higher  energy  neutrons  are  more  damag¬ 
ing,  as  can  be  seen  by  the  D( En )  curve  for  silicon  given  in  References  2  and  5. 

The  "spectral  shape  factor"  <J>  /<J>  given  in  Column  7  shows  the  damage- 

eq 

ability  per  neutron  above  10  keV.  These  values  also  differ  by  nearly  28%, 
in  agreement  with  the  conclusions  of  the  preceding  paragraph. 

Note  that  the  difference  in  <P^/4>  for  exercise  wheel  behind  2"  Pb 
wall  versus  ER-1  and  ER-2  free  field  is  about  1/3  that  for  6"  Pb  +  cave  ver¬ 
sus  free  field,  as  might  be  expected. 

Adding  tissue-shielding  (-3^"  radius)  to  the  6"  Pb  shield  (Run  5) 
results  in  only  a  very  slight  hardening  of  the  spectrum:  <J>eq/$  increases  by 
only  2%.  However,  the  effect  of  the  tissue  on  the  total  neutron  flux  is 
nothing  short  of  dramatic.  While  the  6"  Pb  shield  and  cave  had  no  effect  on 
total  neutron  flux,  the  addition  of  the  3*s"  radius  ohantom  reduced  the  flux 
by  a  factor  of  4!  This  arises  from  the  softening  of  the  neutron  spectrum  by 
lead,  combined  by  the  large  increase  in  hydrogen  cross  section  between  the 
3  MeV  peak  in  E*$(E)  for  free-field  neutrons  and  the  0.8  MeV  peak  for  neutrons 
through  the  6"  Pb  shield.  The  free  field  high-energy  neutrons  are  attenu¬ 
ated  by  about  1%/mm  in  water,  or  a  factor  of  2.6  for  3V'  of  water  (tissue). 

This  is  to  be  compared  to  the  factor  of  4  for  lead-moderated  neutrons  in 
3H"  of  water. 
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Table  5.  Spectral  Parameters  for  AFRRI  TR1GA.  D(En)  Is  the 

total  neutron  "damage"  (energy  deposit)  ir.  Rads  Silicon. 


Note  that  the  measured  total  neutron  flux  above  10  keV,  Table  4, 
and  the  <J>eq,  Table  5,  are  nearly  identical  for  ER-1  and  ER-2  free  field  runs 
(Run  2  and  Run  5,  respectively).  Yet  the  calculated  values  (Table  4)  differ 
by  25%.  This  is  further  established  by  the  independent  nickel-foil  threshold 
neutron  detector  results  shown  in  Column  8  of  Table  5.  This  further  under¬ 
lines  the  necessity  for  normalizing  all  calculational  results  to  the  measured 
data  for  the  measured  configuration  that  is  nearest  to  that  being  calculated. 

The  nickel-foil  monitor  has  proven  to  be  much  more  reliable  than 
the  sulfur  monitor.  For  this  reason,  the  nickel-foil  activation  data  is  pre¬ 
sented  in  Columns  8,  9  and  10  in  Table  5.  Column  8  presents  the  foil  ac¬ 
tivation,  Column  9  the  neutrons  above  10  keV  per  activation/nucleus,  and 
Column  10  the  1-MeV  equivalent  fast-neutron  fluence  per  activation/nucleus. 

The  nickel  foil  is  activated  only  by  neutrons  above  about  3  MeV.  It  is 
interesting  to  compare  the  nickel  foil  activations  for  the  ER-1  free-field 
case  and  6"  Pb  shield  and  cave  case,  where  the  neutron  flux  above  10  keV  is 
the  same.  The  nickel  foil  activations,  on  a  per-kilowatt  basis,  are  a  factor 
of  4  greater  for  the  free-field  case,  which  illustrates  how  much  harder 
this  spectrum  is.  The  biological  dose  and  transistor  damage  for  fast  neu¬ 
trons  are  much  greater  in  the  free-field  case,  where  the  massive  slowing 
down  of  neutrons  has  not  occurred. 

One  can  semi -quantitatively  obtain  neutron-spectrum  "hardness"  data 
of  this  nature  from  as  few  as  two  foils.  The  U-235  data  of  Table  2  show 
nearly  the  same  activation  level  per  kilowatt  for  Runs  1  and  2  (6"  Pb  +  cave 
versus  free  field).  This  agrees  with  the  integral  flux  above  10  keV,  and  it 
shows  almost  no  variation  above  10  keV  whereas  the  nickel  foil  shows  a  fac¬ 
tor  of  4  difference  above  3  MeV.  A  nickel /U-238  "spectrometer"  such  as  this 
is  about  the  simplest  imaginable,  conceptually.  Yet,  it  produces  sensitive 
radiation-hardness  data  for  both  electronics  damage,  and  biological  damage, 
and  possibly  for  damage  to  structural  materials  and  electrical  insulators 
as  well. 
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6.0 


SUMMARY  AND  CONCLUSIONS 


The  radiation  fields  for  the  AFRRI  reactor  have  been  measured  with 
good  accuracy  for  all  energies  of  interest,  and  for  five  widely  different  geo¬ 
metries  commonly  used  at  that  facility.  The  silicon  damage  can  be  determined 
from  these  data  with  an  estimated  precision  of  1-3%.  The  accuracy  has  perhaps 
an  order  of  magnitude  smaller  uncertainty  than  achievable  in  the  1960's,  when 
the  last  determination  was  made.  This  for  the  most  part  is  due  to  increased 
accuracy  in  cross  sections,  to  better  known  nuclear  data  in  other  respects 
(fission  yields,  gartma  fluorescent  yields),  better  calculational  accuracy,  the 
use  of  the  calculated  spectral  shape  in  unfolding  of  the  measured  neutron  spec¬ 
trum,  the  availability  of  large  computers  and  a  neutron-spectrum  unfolding 
code  that  requires  these  machines,  and  the  availability  of  much  better  gaima- 
ray  detectors  and  accurately  calibrated  CNBS)  ganma-ray  line  sources. 

The  calculated  data,  when  normalized  to  the  measured  data  of  nearly 
the  same  configuration,  can  be  of  great  value  for  accurate  estimates  of  neu¬ 
tron  fields.  In  this  way,  accurate  planning  future  experiments  will  be 
greatly  facilitated  if  these  data  are  fully  utilized. 
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